Monte Carlo Study of Maghemite Nanoparticles 
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We consider a simple model of maghemite nanoparticles and study their magnetic properties using 
Monte Carlo methods. The particles have a spherical geometry with diameters ranging from 3 nm to 
8 nm. The interior of the particles consists of core spins with exchange interactions and anisotropy 
given by the values in the bulk material. The outer layer of the particles consists of surface spins 
with weaker exchange interactions but an enhanced anisotropy. The thermal behaviour of the total, 
core and surface magnetizations are calculated as well as the hysteresis loops due to the application 
of an applied field. The effect of the surface anisotropy on the blocking temperature, the coercive 
and exchange bias fields is studied. 

PACS numbers: 75.10.Hk, 75.40.Mg, 75.50.Tt 
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I. INTRODUCTION 

Nanomagnetism focuses on the magnetic behaviour of 
individual building blocks of nanostructured systems as 
well as on combinations of individual building blocks 
that display collective magnetic phenomena. A funda- 
mental understanding of nanomagnetism will lead to the 
development of integrated systems with complex struc- 
tures and architectures that possess new functionalities. 
Proximity effects allow multi-component composites to 
behave as new materials that embrace properties that 
are often mutually exclusive and thus not found in sin- 
gle component systems. Competing interactions and the 
presence of low-lying energetic states and quantum fluc- 
tuations help create the complexity that gives rise to 
unanticipated phenomena in magnetic nanosystems. 

Magnetism and magnetic materials have been tradi- 
tionally studied with phenomenological models. These 
models either work well or must be supplemented by new 
terms in the model to account for unexplained effects. 
Such an approach has limitations at the nanoscale. The 
model may not be able to explain characteristics that de- 
pend on the details of the system at the nanoscale and 
may not have predictive capabilities. Magnetic proper- 
ties at interfaces and surfaces, which make up a large 
fraction of nanostructured and confined materials, are 
quite different from the bulk systems upon which many 
simple models are built. 

Fundamental to understanding the magnetic behaviour 
is the evolution of the magnetism as the structural scale 
descends from the bulk to the nanoscale. Due to reduced 
symmetry, the magnetic anisotropy can be orders of mag- 
nitude larger than in the bulk. This result can lead to 
magnetic frustration and reorientation of the magnetiza- 
tion at the surface and interface. Understanding the com- 
plex atomic spin structure of magnetic nanostructures 
using computational approaches is thus essential to the 
mastering of nanomagnetism itself. Recently, magnetic 
measurements^ on a collection of iron-oxide nanoparti- 
cles arranged on a macroscopic three-dimensional fee ar- 



ray have been observed to exhibit a low coercivity that 
is weakly temperature dependent with no superparam- 
agnetism up to 400 K, whereas only dynamical freez- 
ing and superparamagnetism is observed for a randomly 
packed configuration of the same nanoparticles. A nu- 
merical study^ of the effects of anisotropy, magnetostatic 
interactions and temperature on hysteresis loops for a 
collection of magnetic dipoles on a fee lattice using fi- 
nite temperature micromagnetic simulations found that, 
although general loop shapes resulting from the model 
are in agreement with these experimental results, some 
crucial features are not reproduced, especially as a func- 
tion of temperature. The intrinsic temperature depen- 
dence of the core and surface magnetizations of individ- 
ual nanoparticles needs to be incorporated to correctly 
describe the experimentally observed reversal processes 
in ordered arrays of nanoparticles. In this work we only 
consider the magnetic properties of single nanoparticles 
in order to determine which properties might be impor- 
tant to understand ordered arrays of nanoparticles. 



II. EXPERIMENTAL FACTS 

The magnetic properties of nanoparticles are greatly 
influenced by surface effects due to their finite size^i^ 
and the relative importance of surface sites increases as 
the size of the particle decreases. With decreasing size, 
surface effects lead to a decrease of the ordering tem- 
perature and an enhanced surface anisotropy which dif- 
fers significantly from that of the particle core. Recent 
experiments^ on dilute suspensions of mono-disperse 7- 
Fe203 nanoparticles with an average diameter of 7 nm 
have suggested that the particles block at a temperature 
Tb ~ 20 K which is much lower than that where the core 
develops ferrimagnetic order. Below Tb the particles ex- 
hibit enhanced coercivity and exchange bias. 

With decreasing size, the average magnetic coordina- 
tion number of a 7-Fe203 crystallite is significantly re- 
duced, and finite-size effects become significant. A single- 
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domain ferrimagnetic structure exists in the core, and 
the large number of Fe ions in the unit cell make this 
system sensitive to vacancies distributed throughout the 
octahedral sites and especially at the surface. This leads 
to a magnetic "roughness" which can be thought of as 
an effective surface magnetocrystalline anisotropy that 
is different from the magnetocrystalline anisotropy of 
the nanoparticle single domain core. Hence the single- 
domain ferrimagnetic core of the 7-Fe203 nanoparticle is 
magnetically ordered at a much higher temperature than 
the surface spins. 
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FIG. 1: Saturation magnetization, Ms{T)/Ms{Q), of 7-Fe203 
nanoparticles as a function of temperature^. 

A prominent experimental feature of the 7 nm 7-Fe203 
particle nanomagnetism is the unique temperature de- 
pendence of the saturation magnetization, Ms, shown 
in Fig. [T] The effect of surface spin freezing in the 7- 
Fe203 nanoparticles displayed by Ms{T) in Fig. [l]is well 
described by the jahenomenological description using a 
modified Bloch law 

Ms{T)/Ms{Q) = (1 - A){1 -BT-^ )+ Ae7*~s (1) 

which describes Mg over the complete range of tempera- 
tures. Fits to the data in Fig. [l]with Eqn.([T]) (shown by 
the solid line) identifies a surface spin freezing tempera- 
ture Tf = 15 ± 1 K^. 



III. MODEL 

Bulk maghemite (7-Fe203 ) is derived from magnetite 
(Fe304) which crystallizes in an inverse spinel structure 
with 8 Fe'^^ ions located at tetrahedral sites (T), 16 Fe'^'*' 
ions in octahedral sites (O), and 32 O'^" ions per unit 
cell^"— . The lattice constant of the unit cell is .83 nm. 
Maghemite is formed by the removal of ^ of the Fe atoms 
and vacancies are distributed throughout the O sites to 
ensure charge neutrality. Hence there are 16 — ^ = ^ 
Fe'^+ atoms on the O sites in a unit cell. The ground state 
configuration of the Fe'^"'" ion is which has S = — 
and gj = 2. Competing superexchange interactions 



between the Fe"^"*" T and O sites leads to a ferrimagnetic 
ordering in the bulk material at 1020 K. The moments on 
the T and O sites are anti-parallel and the net moment 

per iron atom is ( io.c, )^I^B = ■'25gjSiiB = 1.25/is. 

In order to study the magnetic properties of maghemite 
nanoparticles, we shall use the following Hamiltonian 

i<Cj i^s i^c i 

where the first term describes the exchange interactions 
with the nearest magnetic neighbours of each iron site, 
the second term describes the single-ion uniaxial sur- 
face anisotropy, the third term accounts for the uni- 
axial core anisotropy and the last term is an applied 
magnetic field. The unit vector is the local normal 
to the surface of the particle. In the spinel structure 
the tetrahedral sites have 12 nearest neighbours on the 
octrahedral sites and 4 nearest neighbours on tetrahe- 
dral sites whereas the octrahedral sites have 6 nearest 
neighbours of each type. The corresponding superex- 
change interactions were taken to be antiferromagnetic 
with the values Jtt = —i2.0K, Jto = — 56.2ii', and 
Joo — —17-2K. These values were taken from the lit- 
erature for bulk maghemitoiir— . Since the electronic, 
structural and magnetic properties of nanoparticles are 
modified in the surface region, we expect an enhanced 
anisotropy and weaker exchange. However, the values 
of these parameters are not well known. For this reason, 
we make the reasonable assumption that the surface-core 
exchange interactions are all weaker by factor of two and 
the surface-surface interactions are weaker by a factor of 
ten. The single-ion site surface anisotropy was given 
the values 1.0,5.0 and 10.0 , while the core anisotropy 
Kc = 0.02. The spins at each site are taken to be classi- 
cal spins of unit magnitude and vacancies on the octrahe- 
dral sites are given a spin magnitude of zero. The oxygen 
ions are considered to be non-magnetic and only serve to 
provide superexchange interaction pathways between the 
occupied Fe sites. The actual spin of each iron atom is 
S" = I and is large enough to be treated classically at fi- 
nite temperatures. However, important quantum effects 
are present at low temperatures. Since our model uses 
classical unit vectors for the Fe ions whereas the actual 
spin length is J — gjS — 5, there is an extra factor of 
J(J -I- l)/3 = 10 in the temperature scale that has not 
been included in the definition of our exchange constants. 

We take the maghemite particles to have a spherical 
geometry with a diameter of L unit cells as shown in 
Fig. [2j The Fe^^ ions in the outermost shaded cells are 
called surface sites and the ions in central cells are core 
sites. For L = 5 the diameter is about 4.2 nm and the 
particle contains 1112 surface sites, 311 core sites and 
170 vacancies. For this size of particle, the majority of 
sites are on the surface. 

We study particles with sizes ranging from L = 4 to 
L = 9 (3 to 8 nm diameter). The number of surface, 
core and vacant sites for L = 4 to 9 are given in Table 
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FIG. 2: Geometry of nanoparticles 
TABLE I: Number of Surface, Core and Vacant Sites 



L Diameterinm) Surface Core Total Vacancies 



4 


3.3 


667 


94 


761 


80 


5 


4.2 


1112 


311 


1423 


170 


6 


5.0 


1683 


751 


2434 


297 


7 


5.8 


2467 


1411 


3878 


473 


8 


6.6 


3358 


2440 


5798 


687 


9 


7.5 


4327 


3891 


8218 


971 



|T1 The vacancies are distributed uniformly throughout 
the octahedral sites in both the core and surface and 
account for approximately 1/9 of the total sites. From 
Table U we see that the majority of the sites arc on the 
surface. For a perfect ferrimagnetic order where the oc- 
cupied tetrahedral sites have their spins anti-parallel to 
those on the octrahedral sites, we can use the data in Ta- 
ble H] to predict the values of the core, surface and total 
magnetization per core, surface and total number of sites 
respectively as shown in Table HIl 

TABLE IL Predicted Core, Surface and Total Magnetizations 



L 


TfLcore 


TTT'Sur face 


mtotal 


4 


.255 


.262 


.261 


5 


.209 


.277 


.262 


6 


.252 


.247 


.248 


7 


.256 


.275 


.268 


8 


.250 


.265 


.259 


9 


.271 


.241 


.255 



IV. MONTE CARLO RESULTS 

We will use Monte Carlo techniques to study the mag- 
netic properties of these nanoparticles. There have been 



many previous theoretical studies of the magnetic proper- 
ties of nanoparticles. The group of Restrepo et. alJ^"— 
have used the Metropolis algorithmic to study the ef- 
fect of surface anisotropy on the sublattice magneti- 
zations with uniform exchange interactions throughout 
the particles. Trohidou et. ali^^— have also used the 
Metropolis algorithm to study nanoparticles with a fer- 
romagnetic core surrounded by an antiferromagnetic sur- 
face layer. Iglesias et. al<^"— have studied both sys- 
tems with a ferromagnetic core surrounded by an anti- 
ferromagnetic surface layer as well as a more realistic 
model of maghemite with the inverse spinel structure. 
Several other groups^S"— have also used the conven- 
tional Metropolis algorithm to study the effects of sur- 
face anisotropy on the magnetic properties of maghemite 
nanoparticles. 

We employ a heat bath algorithm that is a modified 
version of the usual heat bath method used to study 
classical spins to allow for the finite length of the Fe 
spins in maghemite. The modifications only affect the be- 
haviour of the thermodynamic quantities at low temper- 
atures where classical models fail. The heat bath Monte 
Carlo method^i^ generates a sequence of configurations 
which simulate a canonical ensemble of states in thermal 
equilibrium at a constant temperature T. Given a con- 
figuration of the spins in a system, the direction of each 
spin is updated assuming that the spin is in contact with 
a heat bath that puts it into local equilibrium with the 
surrounding spins. At each Monte Carlo step, there is a 
spin update and hence no moves are rejected. This is a 
very efficient way to sample the accessible microstates. 
In the case of classical spins, the spin direction with re- 
spect to the local field H^^'^ can be described using the 
azimuthal and polar angles (f>i and 9i respectively. The 
heat bath algorithm usually allows the cone angle 9 to 
have any value in the range [0,7r]. As shown below in 
Fig. [21 a semi-classical picture of the possible spin ori- 
entations in a field is a set of cones with the local field 
along the cone axis 

We use the heat bath algorithm to generate a value 
of cos 9 but we then place the spin direction relative to 
the effective field onto the nearest cone. This procedure 
has little effect at high temperatures but does modify 
the low temperature dependence of the magnetization 
and energy. We start the simulation at high temperature 
and cool in either a zero or applied field. We measure 
the energy, core and surface magnetizations as a function 
of temperature. At low temperatures we also measure 
hysteresis loops. Typically we use 30, 000 mcs for our 
measurements at each value of H and T. We have studied 
particles with sizes ranging from L = 4 to 9 (3 to 8 nm 
diameter) . 



A. Thermal Dependence of the Magnetization 

The temperature dependence of the magnetization can 
give important information about the properties of the 
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nanoparticles. In the classical model we have a unit vec- 
tor at each Fe site in the nanoparticle. The core, surface 
and total magnetizations are given by 



Total Magnetization 



mtot 



- S 

2^ core 

/ - ^ E 



(3) 



ntot 



ietot 



where ncore,nsurf and ntot are the number of occupied 
core, surface and total sites as given in Table HI We have 
calculated average values of each component of the mag- 
netization as well as the magnitude squared at each tem- 
perature. The particle is initialized at high temperature 
with the spins in random directions. The Monte Carlo 
updates are then performed for 1000 Monte Carlo steps 
(mcs) and the measurements are collected over the next 
5000 mcs at this temperature. The temperature is then 
lowered and the process is repeated. 

The upper panel of Fig. U shows a plot of the magni- 
tude of the total magnetization for particle sizes ranging 
from L = 4toL = 9asa function of T in the zero 
field cooled case. The surface anisotropy constant has 
the value Ks = 5. At high T, the magnetization magni- 
tudes do not vanish but rather approach ^^^^1— with the 
larger values of L lying below those of smaller L. This is 
a finite size effect. At lower temperatures, the magneti- 
zations with smaller L lie below those of larger L. The 
crossover for a given pair of sizes occurs at temperatures 
in the range from T = 50 to T = 70. For the largest 
size pair, the crossover is near Tc ~ 68. We identify this 
temperature with the ordering temperature of the parti- 
cle where strong correlations between neighbouring spin 
directions develop. This value of Tc differs from the ex- 
perimental value of Tc for maghemite particles of these 
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FIG. 4: Magnitudes of the total magnetization (upper panel), 
the core (middle panel) and the surface (lower panel) magneti- 
zations as a function of temperature for particle sizes ranging 
from L — 4 to L — 9. The surface anisotropy constant has 
the value Ks — 5 and the applied field H = 0. 
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FIG. 3: Semiclassical description of spin with S — 5/2 
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sizes by an order of magnitude because our model uses 
classical unit vectors for the Fe ions whereas the actual 
spin length is J — gjS — 5. As mentioned previously, in- 
cluding this factor would yield an ordering temperature 
of the order of 680 K which does lie in the range of the 
experimental values. 

Below Tc, the total magnetization increases with de- 
creasing temperature and displays an even more rapid 
increase below about T ^ 5. The total magnetization in- 
cludes contributions from both the core and surface sites. 
Although the core develops order at Tc, the surface sites 
can remain disordered until much lower temperatures. 
The Monte Carlo approach allows us to separate these 
two contributions. The middle and lower panels of Fig. 
Slshow the magnitudes of both the core and surface mag- 
netizations as a function of T. The core magnetizations 
begin to develop at about Tc ~ 68 and saturate at low T. 
The values of the magnetizations at T = are slightly 
less than those predicted in Table ITll and indicate we have 
an imperfect ferrimagnetic order in the core. In contrast, 
the surface magnetizations remain zero down to very low 
temperatures and only increase substantially only below 
T ^ 5. The particles with small L are dominated by sur- 
face sites and only have a small number of core sites. The 
particles have ferrimagnetic order with a small amount 
of disorder due to the oxygen vacancies on the octrahe- 
dral sublattice which weaken the exchange interactions. 
The surface sites are also affected by the reduced coor- 
dination number at the spherical surface. Note that the 
T = core magnetizations are reduced much more for 
small L and that the surface magnetizations are reduced 
even more than those of the core for all sizes. 

Additional information can be obtained by measuring 
the vector components of the total magnetization as the 
temperature is reduced. Fig. [S] shows the components 
and the magnitude of the total magnetization as a func- 
tion of T for L — 7 at low temperatures. At higher 
temperatures T > 70, all the components average to 
zero which indicates that the particles are paramagnetic 
in this temperature range. Below Tc ~ 68, the com- 
ponents of rhtot fluctuate from one temperature to an- 
other with an increasing amplitude and the magnitude 
increases smoothly. This indicates that the particle is 
developing ferrimagnetic order but it is not blocked in 
this range of temperature. Although the particle has a 
net moment, it is superparamagnetic since the direction 
of the net moment can overcome thermal barriers and 
change its direction spontaneously. At much lower tem- 
peratures (upper panel), the components of rhtot cease to 
fluctuate and become blocked below the blocking temper- 
ature Tb ^ 3. 

The middle and lower panels of Fig. [S] show the compo- 
nents and magnitudes of the surface and core magnetiza- 
tions at low T for L = 7. In both cases, the magnitudes 
vary smoothly but the components fluctuate from one 
temperature to another until below Tb ^ 3 where both 
the surface and core magnetizations become blocked and 
cannot overcome the thermal barriers. The surface mag- 
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FIG. 5: Components and magnitude of the total (upper 
panel), surface (middle panel) and the core (lower panel) mag- 
netization as a function of temperature for particle size L — 7. 
The surface anisotropy constant has the value Ks — 5 and the 
field H = 0. 
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netization appears to develop at a slightly higher tem- 
perature Tg ^ 6 which is about 10% of T^. All of these 
results are for a value of the surface anisotropy Kg = 5. 
Similar behaviour is observed for all particle sizes. 

We have repeated the calculations with a smaller val- 
ues of the surface anisotropy constant. For Kg = 1, the 
components of the magnetizations exhibit strong fluctua- 
tions down to much lower temperatures than for Kg = 5. 
The surface ordering temperature is reduced to Tg ~ 3 
and the blocking temperatures are reduced to Tg ~ 0.5. 
For larger values of the surface anisotropy, the surface 
ordering and blocking temperatures are increased. For 
Kg = 10, the blocking temperatures are increased to 
Tb ~ 4.5 and the surface appears to order at ^ 8. 
This behaviour is again characteristic of all particle sizes. 
For all values of Kg and L, the core and surface magne- 
tizations appear to block at the same temperature but 
have different ordering temperatures. 

TABLE III: Blocking Temperature Tb for various L and Ka 



L Ks = 1 Ks ^5 Ks = 10 



4 


0.3 


2.9 


4.0 


5 


0.4 


3.2 


4.2 


6 


0.4 


3.0 


4.4 


7 


0.5 


3.0 


4.5 


8 


0.5 


2.7 


4.5 


9 


0.5 


3.0 


4.5 



The values of Tb estimated from the temperature de- 
pendence of the magnetization components are tabulated 
in Table IIIII for the different values of L and Kg . The 
blocking temperature seems to be fairly independent of 
L, which is a measure of the particle size, but increases 
with Kg. The blocking temperature of a single domain 
particle with an effective anisotropy K^f f and volume V 
is expected to be proportionaP"*'"^^ to KeffV . The ef- 
fective anisotropy is usually assumed to be both temper- 
ature and size independent and hence one might expect 
Tb to increase with particle size. However, several stud- 
ies on nanoparticles^^""— indicate that Keff increases as 
the particle size decreases. The decrease in Kgff with 
size can compensate for the expected increase with V 
and lead to a size- independent value of Tb- A decrease 



in Keff with V could also be due to stronger temper- 
ature renormalization in smaller particles. Our results 
indicate that the effective energy barrier is independent 
of the particle size over the range of L studied. 

The surface magnetizations for Kg = 5 shown in the 
lower panel of Fig. |4] can be fit to a expression of the form 

(T) = aexp{-T/Tg) + b (4) 

where b describes the finite size effects at high T and for 
each L has the value , ^ The fits were carried out 

over the temperature range [1 , 50] and the fitting param- 
eters are tabulated in Table IIVI The values of Tg are 
slightly larger than the values of Tb determined from the 
components of the magnetization and indicate the tem- 
perature where the surface spins develop order but are 
not yet blocked. 

TABLE IV: Fit to msurf{T) = ae"^/^- + b for A, = 5. 
T a b Ts 

4 0.13 ± 0.01 0.034 ± 0.094 2.91 ± 0.39 

5 0.13 ± 0.01 0.033 ± 0.001 5.83 ± 0.40 

6 0.13 ± 0.01 0.024 ± 0.001 4.75 ± 0.35 

7 0.15 ± 0.01 0.020 ± 0.001 4.69 ± 0.32 

8 0.15 ± 0.01 0.019 ± 0.001 5.93 ± 0.35 

9 0.13 ± 0.01 0.013 ± 0.001 4.29 ± 0.36 

However, if we increase the number of temperature 
points in the range [0.2,10], then the exponential form 
is not appropriate. The modified form 

rUgurfiT) = aexp(-(T/T,)") + 6 (5) 

works extremely well. Table |V] shows the values of Tg 
and the exponent n for Kg = 1,5 and 10. The value of 
Tg increases with Kg in much the same way as Tb but is 
fairly independent of L. Tg is larger than Tb which seems 
to indicate that the surface develops order before the 
particle becomes blocked as the temperature is reduced. 
The exponent n is also independent of L but increases 
with the value of Kg. 





TABLE V: Fit to msurf{T) = ae~'-^^" 


6 for Ks = 1, 5, 10 




Ks = l Ks = 5 


Ks = 10 


L 


Ts n Ts n 


Ts n 



4 2.92 ± 0.01 2.30 ± 0.02 4.40 ± 0.01 3.02 ± 0.04 5.99 ± 0.02 3.75 ± 0.05 

5 3.25 ± 0.03 1.84 ± 0.04 4.82 ± 0.03 2.60 ± 0.05 6.19 ± 0.03 3.32 ± 0.05 

6 3.38 ± 0.02 1.96 ± 0.04 4.99 ± 0.02 2.86 ± 0.04 6.34 ± 0.02 3.72 ± 0.04 

7 3.34 ± 0.02 1.88 ± 0.03 4.88 ± 0.02 2.64 ± 0.03 6.24 ± 0.01 3.27 ± 0.02 

8 3.64 ± 0.02 1.81 ± 0.03 5.19 ± 0.02 2.54 ± 0.03 6.57 ± 0.01 3.19 ± 0.02 

9 3.59 ± 0.02 1.92 ± 0.03 5.15 ± 0.01 2.69 ± 0.02 6.54 ± 0.01 3.27 ± 0.02 
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Both Tb and Ts display the same quahtative behaviour 
with Ts > Tb as a function of L and Kg- This suggests 
that the particle blocking is associated with the surface 
ordering and that the surface anisotropy plays an impor- 
tant role in both quantities. 



B. Field Cooling 

We have also studied the magnetic properties of the 
nanoparticles when they are cooled in the presence of an 
applied field. The field cooling procedure starts from a 
random spin configuration at high T and then the tem- 
perature is slowly reduced to T = 1. During the cooling, 
we measure the total, core and surface magnetizations. 
Fig. [6] shows the magnetizations for both the zero field 
and non-zero field cases with L = 7 and Kg — 5. The 
lower panel of the figure shows that the magnetizations 
are increased in the field cooled case. At the lowest tem- 
perature, the surface magnetization has almost doubled 
compared to the zero field case. The magnitude of the 
applied field is the same as the magnitude of the surface 
anisotropy and these contributions to the energy of the 
particle are in competition. 

The surface anisotropy favours the surface moments to 
be normal to the surface whereas the applied field favours 
alignment in a particular direction. When cooled in zero 
field, the surface anisotropy leads to some of the surface 
moments pointing outward from the surface and some in- 
ward as both satisfy the normal arrangement. Hence the 
surface anisotropy can lead to a net moment reduction. 
In the presence of a field, the direction of the surface 
moments will change and will point outward or inward 
according to whether the outward normal has a compo- 
nent parallel or anti-parallel to the field. The net effect is 
an increase in the net surface moment as shown schemat- 
ically in Fig. [71 



Magnetizations for L=7 and H=0 
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FIG. 6: Magnitudes of the total, core and surface magneti- 
zations as a function of temperature for particle size L — 7 
when cooled in a field of magnitude H = (upper panel) and 
a field H = 5 (lower panel). The surface anisotropy constant 
has the value Kg = 5 in both cases. 



C. Hysteresis Loops 

After cooling the particle in either zero field or non- 
zero field, we have measured the hysteresis properties of 
a nanoparticle. In each case we begin at the lowest tem- 
perature with a field applied in the ^-direction and slowly 
reduce its magnitude to zero and then slowly increase the 
field in the reversed direction until the magnitude is the 
same as it was at the beginning. This process is then 
carried out in the reverse direction. At each value of the 
applied field, we perform 30,000 mcs to measure the com- 
ponents of the magnetization and the energy. Magnetic 
contributions to the total magnetization per site arising 
from the surface and core sites were computed separately. 

Fig. [S] shows typical hysteresis loops for L = 7 and 
Ks = 5 for the zero-field cooled case at various tem- 
peratures. The upper panel corresponds to the range 
1 < T < 5 and indicates that the width of the hystere- 
sis loop quickly narrows as the blocking temperature is 



reached. Tb was estimated to be ~ 3 for this size of 
particle. The middle panel shows hysteresis loops in the 
narrower range 1 < T < 1.4. As the applied field is re- 
duced in magnitude and reversed, the component of the 
magnetization in the z-direction slowly decreases until 
it rapidly changes direction near ^ —0.2. However, 
the reversal is not complete until ^ —0.4. Similar 
behaviour is observed along the return path as is in- 
creased to positive values. At the lowest temperature 
T = 1, the coercive field (as measured by when the mag- 
netization changes sign) is He ~ 0.25. The T = 1 loop 
appears to have sub-loops when the magnitude of the 
applied field is between 0.2 and 0.4 which indicate that 
complete reversal does not occur immediately. 

The lower panel in Fig. [5] shows the energy per particle 
as a function of the applied field for T = 1, 1.2, 1.4. As 
the field is reduced from — 0.6, the energy increases 
reaching a maximum near Hz ~ — 0.2 and drops abruptly 
as the field is reduced further. A second abrupt decrease 
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FIG. 7: Zero-field cooled and non-zero field cooled configura- 
tions 



occurs near Hz ^ —0.4. These abrupt decreases indi- 
cate the hmits of metastable configurations of the parti- 
cle as the applied field is changed. As the temperature 
increases, the size of the abrupt energy jumps decrease 
and their location moves towards Hz = 0. Hence the 
reversal appears to be a two step process. 

Fig. [9] shows the hysteresis loops for a particle which 
has been cooled in a field H — 5 applied in the z- 
direction. The behaviour is very similar to that in the 
zero-field cooled case shown in Fig. [8] except that the 
hysteresis loop appears to have a small shift to the left 
and is no longer symmetric about Hz = 0. This feature 
indicates that a small amount of exchange bias may be 
present. However, this shift disappears as soon as the 
temperature is increased. Both the core and surface ex- 
hibit the two stage reversal process for this value of Kg- 
This suggests that the particle reversal is quite different 
from a single domain reversal mechanism where the core 
reversal occurs separately and is affected by the surface. 
On the contrary, different regions of the particle have 
their core and surface acting together at two different 
reversed fields. 

We have also studied the hysteresis properties for 
smaller and larger values of the surface anisotropy con- 
stant Kg. Fig. [TOl shows results for L = 7 and Kg = 10 
after cooling in a field. The upper panel shows ruz at sev- 
eral low temperatures where a large coercivity is present. 
The two step reversal process found for = 5 now seems 
to be replaced by a continuous process which occurs over 
a wide range of reversal fields. The lower panel shows the 
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FIG. 8: Hysteresis loop for L — 7 and Kg = 5 in the zero 
field cooled case. The upper panel shows the z-component of 
the total magnetization as a function of the applied field in 
the range 1 < T < 5. The middle panel shows the behaviour 
over a narrower temperature range. The lower panel shows 
the energy as a function of field for T = 1, 1.2, 1.4. 
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Hysteresis Loop for L=7 H=5 and Kg=5 
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FIG. 9: Hysteresis loop for L = 7 and Ks = 5 in the non-zero 
field cooled case. 



hysteresis loop for both the core and surface magnetiza- 
tions when Kg = 10 at the lowest temperature T = 1. 
The core and surface magnetizations reverse at different 
fields. As the applied field is reduced from it's maximum 
value, the surface magnetization decreases linearly un- 
til the point where the core magnetization changes sign. 
This is followed by another linear region which has a 
much larger slope. Hence, for larger values of the surface 
anisotropy, the reversal of the particle's magnetization 
again appears to be a two stage process but it is quite 
different from that at smaller values of Kg- 

Fig. [TT] illustrates schematically the two stage reversal 
process for the case Ks = 5. In frame (a), the core and 
surface magnetizations are polarized in the direction of 
the applied field. As the field is reduced and reversed 
as in frame (b) , the magnetizations of both the core and 
surface in the equatorial region of the particle develop 
components perpendicular to the field. Note that the 
figure is a cut through the particle in the x — z plane. 
Since the anisotropy favors the moments to be normal 
to the surface, it is this region which can overcome the 
applied field first. As the field is reversed further as in 
frame (c), the core and surface spins in the polar re- 
gions reverse but the equatorial regions maintain a large 
transverse component. Finally in frame (d), the core and 
surface magnetizations are polarized in the reverse direc- 
tion. Similar reversal processes have been reported by 
Berger et. al»^ 

Fig. [T^ shows a similar picture for the case when 
Kg = 10. Stages (a) and (d) are the same as in the 
previous case but stages (b) and (c) are different. In (c), 
the core magnetization reverses but the surface magne- 
tization does not. The surface magnetization reverses 
when the applied field is reduced further as in (d). 



Hysteresis Loop for L=7 H=5 and K3=10 
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FIG. 10: Hysteresis loops for L = 7 and Ks — 10 at low 
temperatures (upper panel) in the field cooled case and the 
core and surface contributions at T = 1 (lower panel), at 
T = 1. 



V. DISCUSSION AND CONCLUSION 

In this paper we have used Monte Carlo methods to 
study the magnetic properties of maghemite nanoparti- 
cles. The Monte Carlo method employed is a modified 
version of the usual heat bath method used to study clas- 
sical spins to allow for the finite length of the Fe spins in 
maghemite. The modifications only affect the behaviour 
of the thermodynamic quantities at low temperatures 
where classical models fail. At high temperatures the 
core spins begin to develop a ferrimagnetic order with a 
net moment. However, the direction of the moment is 
fluctuating and the particle exhibits superparamagnetic 
behaviour down to much lower temperatures. The sur- 
face spins do not develop ferrimagnetic order until a much 
lower temperature ^ 5. The net moment of the parti- 
cle continues to fluctuate until the blocking temperature 
Tb ^ 3 is reached and a coercivity develops. The temper- 
ature dependence of the surface magnetization obtained 
from our Monte Carlo results is very similar to that ob- 
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FIG. 11: Schematic description of the magnetization reversal process for Ks — 5. 




h| la) H| Hj "1 

FIG. 12: Schematic description of the magnetization reversal process for = 10. 



served experimentally^. 

Both the surface ordering temperature and the block- 
ing temperature do not depend on the size of the 
nanoparticle. However, they both increase approxi- 
mately linearly with the surface anisotropy constant K^- 
The fact that Tb does not depend on particle size in 
the range of sizes studied indicates that the effective 
anisotropy barrier is temperature dependent and leads 
to an effective anisotropy K^ff ~ 1/^- Increases of 
Kf^ff with decreasing nanoparticle size have been ob- 
served experimentally^"—. 

When the particles are cooled in an applied field, the 
magnetizations increase. The largest change occurs for 
the surface magnetization. Hysteresis loops were studied 
in both the zero field and non-zero field cooled cases. 
The differences in the loops for these different cooling 
histories were quite small. Experimental systems exhibit 
a shift of the hysteresis loop when cooled in a field. This 
shift corresponds to an exchange bias field. We find little 
evidence of exchange bias in our model except at the 
lowest temperature. 

The surface anisotropy has a strong effect on the hys- 
teresis loops and the nature of the moment reversal of 
the particles. For small values of Kg, the blocking tem- 
perature is reduced significantly and the hysteresis loops 
only display superparamagnetic behaviour with no coer- 
civity. The particles reverse their magnetization by fol- 
lowing the field direction. For intermediate values of Ks, 
the blocking temperature is increased and the hysteresis 



loops exhibit coercivity. There is a main loop as well 
as minor loops. The reversal of the moment appears to 
be a two stage process. At large values of the surface 
anisotropy, the blocking temperature is increased again 
and the hysteresis loops exhibit a ten-fold increase in the 
coercivity. The minor loops are merged with the main 
loop and the moment reversal occurs in two steps with 
the core spins reversing first followed by the surface spins. 

The results presented here^ are our first attempt at 
studying the magnetic properties of nanoparticles. Some 
of the results agree qualitatively with recent experiments 
but others do not. The absence of exchange bias in our 
model could be due to the fact that our model is too 
simple or it could be due to our Monte Carlo algorithm. 
The Monte Carlo algorithm is a method used to study 
equilibrium properties of magnetic models, whereas, hys- 
teresis and exchange bias are determined by metastable 
configurations. Further work should include modifying 
the model or the algorithm or both. One of the prob- 
lems with our calculations is that they require a great 
deal of cpu time. Each hysteresis loop uses 30,000 mcs 
for each of 200 field values at a single temperature which 
typically requires 24 hours to complete. Parallelization 
of the code could speed up this process significantly and 
allow us to vary both the exchange parameters and ap- 
plied field strengths over a wider range. A study of the 
same nanoparticles arranged on a lattice would require 
such improvements. 
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